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Abstract 

Structural and electronic properties of the a- and y-phases of cerium sesquisulfide, Ce2S3, 
are examined by first-principles calculations using the GGA+U extension of density func- 
tional theory. The strongly con^elated /-electrons of Ce are described by a Hubbard-type 
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on-site Coulomb repulsion parameter. A single parameter of U' = 4 eV yields excellent 
results for crystal structures, band gaps, and thermodynamic stability for both Ce2S3 al- 
lotropes. This approach gives insights in the difference in color of brownish-black a-Ce2S3 
and dark red Y-Ce2S3. The calculations predict that both Ce2S3 modifications are insulators 
with optical gaps of 0.8 eV (a-phase) and 1.8 eV (y-phase). The optical gaps are deter- 
mined by direct electronic excitations at k = F from localized and occupied Ce 4/-orbitals 
into empty Ce SJ-states. The /-states are situated between the valence and conduction 
bands. The difference of 1 eV between the optical gaps of the two Ce2S3 modifications 
is explained by different coordinations of the cerium cations by sulfur anions. For both 
Ce2S3 modifications the calculations yield an effective local magnetic moment of 2.6 /jb 
per cerium cation, which is in agreement with measurements. The electronic energy of the 
a-phase is computed to be 6 kJ mol^' lower than that of the y-phase, which is consistent 
with the thermodynamic stability of the two alio tropes. 
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1 Introduction 



Rare earth (RE) elements are essential ingredients in a large variety of materials of 
vital technological, economic, and ecological importance. This includes 
energy-efficient lighting systems, plasma displays, medical imaging, automotive 
catalysts, permanent magnets, and pigments l^. These applications depend on the 
unique structural, electronic, optical, magnetic, and chemical properties of RE 
compounds. The origin for many of these properties is the presence of localized 
4/-electrons interacting with the valence electrons. In this context, the present 
work focuses on cerium sesquisulfide (Ce2S3), which is used as an 
environmentally responsible pigment with high temperature stability [|20. Beyond 
this specific application, Ce2S3 can be seen as a representative of a large class of 
RE compounds containing 4/-electrons. Therefore, the approach and results of the 
present work has implications for many RE-containing compounds and their 
applications. 

The localization of the 4/-electrons involves strong correlation effects which are 
not captured by current standard methods within density functional theory (DFT). 
The result is that these DFT calculations fail to reproduce, e.g. the non-metallic 
ground state of Ce^Ss \%\4] and overestimate the mass density of crystal 

nun 

structures of other cerium compounds I15|,|6|,|7D. Hence, approaches beyond 
conventional DFT methods are required to describe Ce2S3 realistically from first 
principles. The present work demonstrates the usefulness of a DFT method 
extended by a Hubbard-type on-site Coulomb repulsion parameter for treating the 
Ce 4/-electrons (the GGA-i-U approach). 

Cerium sesquisulfide exists in three allotropic forms m,^. The a-phase is the 
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low-temperature modification and has a brownish-black (rusty) hue nlCm . Above 
~900°C a-Ce2S3 transforms irreversibly into the dark red (maroon) colored 
P-phase. The (3-modification has been described as an oxysulfide with the formula 

a. 



CeioSis-jfOx lllin . The oxygen can be partialy rep. 



aced by sulfur. However, a 



complete substitution appears to be impossible Illin . despite other assertions of 
Picon et al. jioj . The structure of the |3-phase has been described as complex, but 
has not been resolved. For these reasons, the p-phase is not considered in the 
present computational work. The high-temperature allotrope, Y-Ce2S3, is 



preferentially formed between 1 100°C and 



stabilized at temperatures below 800°C 11121 . 



200°C. However, Y-Ce2S3 can be 



at~2100°C lie, 



IJ. 



The melting point of the y-phase is 



The y-phase of Ce2S3 is a nontoxic inorganic high performance pigment of dark 
red (burgundy) hue, which is used on an industrial scale in the coloration of 



1511 suggest that the dark red 



plastics, paints, and coatings |12|]. Optical studies lll4l 
color of y-Ce2S3 arise from dipole- allowed excitations of electrons from 4/-states 
into the empty 5^^-states (4/^5J° — > 4/°5J^) of the same cerium cation. The three 
modifications of Ce2S3 have an insulating nature, as is already clear from their 
color. 



The magnetic susceptibility of a— Ce2S3 and y-Ce2S3 at room temperature are 
similar u^: 2.9 x 10^ (a-phase) and 2.7 x 10^ (y-phase), in SI units. The y-phase 



exhibits ferromagnetism below 3 K and the measured effective magnetic moment 



17 



1811 . The magnetic moment is 



of each cerium cation is peff = 2.5 
equivalent to that of an isolated Ce^+ ion which has the electron configuration 
[Xe]4/' and the ground state ^^5/2 lll9ll . Thorough literature searches did not 
yield anything about a magnetic ordering of a- Ce2S3 at low temperatures. The a 



and y- phases of Gd2S3 have the same magnetic ordering at low temperatures t2m- 
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Moreover, a-GdiSs and y- Gd2S3 are isomorphic to a- Ce2S3 and Y-Ce2S3, 
respectively. For these reasons, we shall assume in the following that both a- 
Ce2S3 and Y-Ce2S3 are ferromagnetic at low temperatures. 



The conduction band (cb) of y-Cej^i is mainly formed by empty Ce 5 J-states 



w 



lile the valence band (vb) consists predominantly of sulfur 3jc-orbitals 
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2 in . Due to their rather extended character and large energy 



dispersion Ce 6^'- states contribute to both the vb and the cb. The measured gap 
between the upper edge of the vb and the lower edge of the cb of y-Ce2S3, Eg, is 



between 2.7 eV U5 



23|] and 3.3 eV jlSl] . In agreement with group-theoretical 



arguments ll24ll the smaller band gap has been assigned to phonon-assisted 
(indirect) electronic excitations. The larger band gap is associated with direct 
electronic transitions. The temperature dependence of Eg between 250 and 350 K 
is approximately —0.1 /jeV K^^, which has been derived from y-La2S3 and 
Y-Dy2S3 \2^. Both compounds have the same structural topology as y-Ce2S3. The 
corresponding values of typical semiconductors like CdS and CdSe are around —4 
A/eV K 1. 



Photoconductivity spectra ll25ll as well as measurements of the photoluminescence 
I26 I of Y-Ce2S3 indicate that the occupied 4/-states of cerium cations are 
energetically located between the vb and the cb. The gap between the 4/- states 
and the cb obtained from these measurements is in agreement with the optical gap 



15|] . Up to now. 



of 1.8-2.0 eV derived from the dielectric function [|14l 
photoemission spectroscopy studies of Y-Ce2S3 have not been able to resolve 
prec isely the energetic position of the occupied Ce 4/-states with respect to the vb 



UK 



15h . However, the measurements of Mauricot et al. ll25h reveal that the 



occupied Ce 4/-states are energetically 0.8 eV above the vb. 
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Earlier band structure calculations within DFT in the local density approximation 
(DFT-LDA) on model structures for Y-Ce2S3 demonstrated qualitatively the 
importance of 4/-electrons for the optical properties M>\4\- However, these 
calculations have several limitations, which hamper quantitative predictions. These 
limitations mostly stem from the well known inability of DFT-LDA methods to 
describe correctly the localized 4/-states. First, when the 4/- states are treated in 
the valence, DFT-LDA calculations predict a metallic (Kohn-Sham) band structure 
for Ce2S3, with /-states at the Fermi level, in contradiction with the insulating 
nature of this material. The fact that, within DFT-LDA, these /-bands are found to 
be split off from the sulfur p-bands is indeed suggestive that the optical transitions 
from 4/ to 5d bands are playing an important role in the optical properties. But it 
is clearly not possible to go beyond this qualitative observation in a pure 
DFT-LDA framework, let alone to reach a quantitative description of the optical 
gaps of Ce2S3 derivatives Second, the models used to describe Y-Ce2S3 in 

Refs. [|3i|4|] did not capture the electrostatics of the real compound adequately. The 
charges of the cations in the model significantly deviate from that of trivalent 
cerium. The deviation was one of the factors for the DFT-LDA calculations to 
underestimate the energy difference between the 4/-states and the vb. The result 
was an overestimation of the overlap between the 4/-states and the vb and an 
artificially large splitting of band states on top of the vb yO. Third, the previous 
LDA-DFT calculations did not perform structural relaxation studies or attempt to 
predict crystal structures. However, the electronic structure and, hence, optical 
properties are very sensitive to interatomic distances. Furthermore, it is widely 
documented that DFT-LDA and even the generalized gradient approximation 
(GGA), an extensions of DFT-LDA, severely overestimate the mass density of 



/-electron compounds when the /-states are treated in the valence 



M- Our 



own experience with alkali-doped cerium sesquisulfides is that GGA calculations 
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underestimate the unit cell volume of these compounds by as much as 6% BOtl . 



The main goal of the present study is to gain a deeper understanding of the 
relations between the chemical composition, the crystal structure, and the optical 
properties of both a-Ce2S3 and Y-Ce2S3. Our study addresses both the crystal 
structure as well as the electronic properties quantitatively within a single 
computational framework, namely the GGA+U approach. Analogue to the 
LDA+U method 1(1. the GGA+U approach uses the DFT framework but extends 
the GGA description of exchange-correlation effects in the localized 4/-shell by a 
many-body interaction term. This term is then treated in a mean-field like manner, 
allowing for a simple description of strong correlation effects within the Ce 
4/-states. The GGA-i-U approach allows for a full optimization of the 
crystallographic parameters of Ce2S3, and captures the role of the 4/-electrons in 
the optical transitions. As we shall see, Ce2S3 is found to be an /-electron Mott 
insulator, in which the Ce 4/- states are split into lower occupied and upper 
unoccupied sub-bands by the strong on-site Coulomb repulsion. Optical 
transitions responsible for the color of these compounds are between the lower 
unoccupied Ce 4/-states and the bottom of the Ce SJ-like cb. 



2 Crystal structures 



The a-phase of Ce2S3 crystallizes in an orthorhombic lattice (crystallographic 
space group Pnma). The unit cell contains four formula units with two and three 
symmetry inequivalent cerium and sulfur ions, respectively lillll . Figure [IJa) 
illustrates two different kinds of interlocked Ce-S skeletons which are the 
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coordination polyhedra of the four Cel and the four Ce2 cations. The coordination 
polyhedra (Cel)S8 and (Ce2)S7 have the shapes of a bi-capped trigonal prism and 
a rather distorted mono-capped octahedron, respectively. The comers of the 
trigonal prism are S2, S3, and S3' whereas the vertices of the two caps are the 
sulfur ions SI'. The cap tip of the octahedron is the sulfur ion S2'. Two of the 
three symmetry inequivalent sulfur anions (S2 and S3) are bound to five cerium 
cations. The cations form a rectangular pyramid. The five cerium cations 
encompassing the remaining symmetry inequivalent sulfur anion (SI) form an 
axial distorted trigonal bipyramid. The Ce-S distances of the two coordination 
polyhedra are in the range of 2.85-3.14 A with average distances of 3.01 A 
[(Cel)S8] and 2.93 A [(Ce2)S7]. 



The crystal structure of Y-Ce2S3 has a body-centered cubic lattice. It is isomorphic 
to the ThsP4 structure and contains vacancies on the cationic sublattice 



32L 



3311 . The cerium cations are located on the 12a positions which have the 
point group symmetry ^4; the sulfur anions occupy sites with C3 symmetry (16c 
positions). Each cerium cation is coordinated by eight sulfur anions forming a 
triangular dodecahedron as depicted in Fig.[TJb), with two different 
nearest-neighbor distances [bj: r(Ce-S)=2.89 A and r(Ce-S')=3.09 A. Each 
sulfur anion is located on a threefold axis and surrounded by six nearest cerium 
cations, provided all 12a positions are occupied. The six cations form a distorted 
octahedron. However, when all the 16c and 12a position are occupied with cations 
and anions, respectively, the stoichiometry would be Ce3S4. In order to achieve 
the stoichiometry of y-Ce2S3 one out of nine cations has to be removed and 



vacancies on the cerium sublattice, Vce, have to be created 



3 X [Ce3S4] = CogSn ^ CegVceSiz = 4 x [Ce2S3] 
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Note that both, 'Y:Ce2S3 and 06384 have very similar lattice parameters differing 
only by 0.06% [|32l. l34|l . The tiny volume difference is evidence for a rather rigid 
Ce-S framework. 



Raman spectra of several RE sesquisulfides in the y-phase give evidence for 
randomly distributed vacancies 113 511 . However, simple electrostatic energy 
calculations predicted that a completely random distribution of the voids on the 
cerium sublattice yields an energetically less stable ground state than a vacancy 
ordering such that the resulting crystal structure of y-CeaSs has the symmetry of 
the space group 742 J 113 311 . A highly symmetric ordering of the vacancies is likely 
in well-annealed y-phases whereas rapidly quenched samples might have a 
random distribution of vacancies ll33ll . Vacancies in the sulfur sublattice are rather 
unlikely 



Figure [TJc) depicts the smallest conventional crystallographic supercell that 
accommodates the stoichiometry of y-Ce2S3. The supercell is three times as big as 
the conventional unit cell of Ce3S4 and has the composition 4x [Ce8VceSi2]. The 
distances of a cerium cation to its first, second and third nearest cations in y-Ce2S3 
are 7.5, 4.8 and 4.0 A B4i1 . In this structure the voids on the cation sublattice create 
three symmetry inequivalent cerium cations. The supercell can be transformed 
into a primitive unit cell half as large as the conventional cell shown in Fig. Ulc). 



3 Computational Method 



All GGA+U calculations have been performed with the Vienna ab initio 



simulation package (VASP) [137 



38L139L 



4111 as integrated in the MedeA 



environment [|42i1 . VASP uses the projector-augmented wave (PAW) method [|430 
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39I]. 



in the implementation of Kresse and Joubert [39']. Exchange- and correlation 
effects are described within the GGA by the density functional of Perdew, Burke 



and Ernzerhof 14411 . The PAW potentials take into account 12 valence electrons of 
each Ce {5s^5p^6s^5d^4f^) and six valence electrons of each S (3s^3p^) atom. 
The plane wave energy cutoff has been set to 333 eV and 299 eV for a-Ce2S3 and 
y-Ce2S3, respectively. 



Within the present GGA+U approach 11401 l45h the strong electron correlation 



among the 4/-states is described with the parameter U' — U —J.ln the latter term, 
U is the spherically averaged screened Coulomb repulsion energy required for 
adding an extra electron into the Ce 4/- states and / represents the corresponding 
screened exchange energy. The GGA+U one-electron potential is defined as 
= dEoGA/dpfj + U'{0.58ij - pfj), whereby Egga is the GGA electronic 
enenergy, is the density matrix of /-electrons with spin a , and 5/j is the delta 
function. The double counting correction applied is equivalent to that of the 
fully-localized limit. The case U' = is the limit of pure GGA. 

The parameter U' can be calculated from first principles without knowledge of 



experimental data 114611 . In this work however, instead of performing such 
calculations (which are computationally heavy because of the large numbers of 
atoms in the unit cell), we have investigated whether U' can be chosen in such a 
way to reproduce with reasonable accuracy both the experimentally measured unit 
cell volume, as well as the optical gap of the Y-Ce2S3 allotrope. We found that the 
value [/' = 4 eV achieves this goal. This value is in agreement with Coulomb 



repulsion parameters calculated for other cerium compounds [|47 



The Kohn-Sham equations are solved self-consistently in their spin-polarized form 
employing an iterative matrix diagonalization (blocked Davidson scheme) and the 
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.[0]. 



spin interpolation of Vosko et al. [49]. The convergence criterion to obtain the 
energy is 10^ eV. The integration over the first Brillouin-zone is performed using 



symmetry adapted Monkhorst-Pack grids 15011 of k-point sets of various size. To 
calculate energies, energy gradients (atomic forces as well as the stress tensors), 
and charge densities the k-point sets are 3 x 7 x 3 (a-Ce2S3) and 3x3x3 
(y-Ce2S3). All k-point sets as well as the plane wave energy cutoff has been 
determined within a convergence test to yield to a numerical precision of 5 x 10^^ 
eV. 

The crystal structures of a-Ce2S3 and y-CeiSs as described in Sec. [2lhave been 
optimized such that atomic positions and the crystallographic cell parameters are 
relaxed simultaneously. The optimization procedure has two steps. In the first step, 
the conjugate-gradient algorithm is used and the optimization stops when all 
energy gradients are smaller than 0.02 eV A^^ Then the crystal structures are 
further refined within a quasi-Newton algorithm with a very tight convergence 
criterion of 0.002 eV A^i. 

Energy gradients have been calculated with a Gaussian integration method and a 
smearing of o = 10 eV. The linear tetrahedron method with Blochl corrections 
has been used to obtain accurate energies and the self-consistent charge densities. 
(This integration scheme is inappropriate for calculating energy gradients.) The 
self-consistent charge densities are used to calculate the band structures and the 
densities of states (DOS's). All DOS's have been calculated using the linear 
tetrahedron method without corrections and the same k-point sets applied to 
calculate the energies (see above). 

Since the calculations correspond to zero temperature, the two Ce2S3 phases adopt 
the ferromagnetic state. In all calculations the ferromagnetic state is described 
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with all /-electron spins alligned parallel (spin-up). The band structures and 
DOS's have been created from the charge density of the electrons with spin-up. 
This is sufficient since in the present approach only the spin-up states contain 
occupied / states while the S 3p-vb and Ce 5(i/6i'-conduction states are similar 
for both spin orientations. 

4 Results and discussion 

4. 1 a-Phase 

Table [H summarizes calculated lattice parameters and interionic distances of the 
two interlocked (Cel)S8 and (Ce2)S7 skeletons shown in Fig.[TJa). As announced 
above, the chosen value U' = 4 eV yields a relative error on the unit cell volume 
which is smaller than 1% in comparison to experiments. A comparison with the 
corresponding experimental distances reveals remarkable agreement. Most of the 
calculated distances are less than 0.9% larger than those determined from x-ray 
diffraction measurements. Exceptions are the calculated cell length c and the bond 
lengths r(Cel-Sr) and r(Ce2-S3) which are nearly equal or smaller, respectively, 
than the corresponding measured distances. The discrepancies are almost zero for 
c, -0.1% (r(Ce2-S3), and around -1.2% for the longest of the two Cel-Sl' 
distances. The calculations tend to strengthen one of the Cel-Sl' bonds at the 
expense of the other thereby slightly decreasing the sulfur coordination of Cel. 
Note that the measurements and the calculations yield for the two distances 
r(Cel-Sr) average values of 3.12 and 3.11 A, respectively. The average of the 
calculated distances between cerium and sulfur within the (Cel)S8 and (Ce2)S7 
skeletons are 3.02 and 2.94 A, respectively. Both distances are only between 0.4 
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and 0.7% larger than the corresponding distances of the x-ray structure. Apart 
from r(Cel-Cel) all other calculated Ce-Ce distances are larger than the 
measured values with deviations not larger than 0.5%. 

The energy of formation, Eform, is a measure for the stability and is the difference 
between the energy of the entire system and the energies of the electronic ground 
state of the constitutive elements in their standard state. For a-Ce2S3 the 
calculations yield Eform = —820 kJ mol~^ CeiSs. 

Figure [2] shows the band structure and the corresponding partial DOS (PDOS) 
calculated for the fully optimized crystal structure of a-Ce2S3. The electronic 
bands are drawn along the standard pathway within the irreducible wedge of the 
Brillioun zone of simple orthorhombic lattices. The PDOS illustrates that the vb, 
i.e. the occupied and evenly dispersed electron bands below approximately —0.8 
eV, are mainly formed by S 3p-orbitals with contributions of Ce 3J-orbitals. This 
indicates a significant covalent contribution to the bonding in a-Ce2S3. 

On top of the vb, i.e. the occupied electron bands between —0.8 eV and zero 
energy are two kinds of occupied 4/-orbitals namely those of the cations Cel and 
Ce2, respectively. In the following occupied 4/-states are termed as lower 
Hubbard bands (LHB's). Both LHB's have rather large admixtures of S 3/?- states. 
This indicates a participation of Ce 4/-orbitals in the covalent bonding. The 
PDOS of Fig. |2b) shows that the LHB of the Cel cations is located at a lower 
energy than the LHB of the Ce2 cations. The main reason for the different 
positions must be the different sulfur environments of the Cel and Ce2 cations as 
depicted in Fig. [Ha). The calculations yield almost one /-electron localized on 
each cerium cation. The computed occupation numbers of the 4/-states, n f, of the 
two symmetry inequivalent cerium cations are 1.12 (Cel) and 1.14 (Ce2). 



13 



Occupation numbers of rif > I axe consistent with admixtures of S 3jc-states to Ce 
4/-orbitals. Assuming that the total orbital moment of the electronic ground state 
of all cerium cations is 7 = 5/2 the effective local magnetic moment is around 2.6 
/Js per cerium cation. This is consistent with experiment. 

The LHB's have different features within the chosen part of the Brillouin zone. On 
the one hand the LHB's have significant dispersion (on the scale of ~ 0.5 eV) 
between k = T and the center of the line between k = F and k = X (E line). On 
the other hand the LHB's are rather flat between the center of the E line and k = 7 
and along the line between k = 7? and k = T. The LHB's are almost degenerate at 
k = 7? and k = S. 

The smallest gap between the upper edge of the LHB of Ce2 and the lower edge of 
the unoccupied states (cb) is 0.8 eV and is at k = P. At the same k-point is the 
smallest separation between the vb and the cb (1.34 eV). The cb has 
predominantly Ce 5c^-character with small contributions of S 3p-orbitals. Like the 
LHB's but more distinct, the cb exhibits a strong dispersion between k = T and 
k = X. The PDOS of the cb at its lower edge is tiny and increases in small steps 
between 0.7 and 1.8 eV. 

The inset of Fig.|2b) depicts the presence of Ce 6s- and Ce 6p-orbitals in the 
relevant electron bands of Fig. Oa). Their contributions within the relevant energy 
range around the Fermi energy are relatively small because the 65- and 6p-states 
are spatially extended and energetically dispersed. Inside the cb between 2.0 and 
4.2 eV are unoccupied Ce 4/-orbitals which form the upper Hubbard band (UHB). 

The band dispersions of the LHB and the cb are examined in terms of an 
orbital-projected band-structure (sometimes called "fat-band" representation). 
Figure [3]illustrates contributions of atomic states to relevant energy bands of Fig. 
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l2ta) between k = T and k = X at different k-points. The three lowest unoccupied 
states at k = r are due to anti-bonding interactions between cerium 5d- and sulfur 
3j!7-states. The contribution of S 3p-orbitals to the three states is between 5 and 
15% whereby only states of the atoms S2 and S3 are involved. Admixtures of 
3/>-orbitals of SI are negligibly small. The corresponding bonding states are 
located between —1.9 and —2.3 eV. The splitting between bonding and 
antibonding states of Cel and S3 is between 0.7 and 1.4 eV smaller than that of 
Ce2 and S2. The smaller splitting and, hence, the distinct dispersion of the two 
lowest unoccupied bands is due to smaller overlap between 5d- and 3p-orbitals of 
Cel and S3, respectively. 

The dispersions at the upper edge of the LHB are also due to admixtures of S 
3j!7-orbitals to the Ce 4/- states. The two highest occupied electronic states of the 
LHB contain at k = F both 44% of S 3p-orbitals whereby all sulfur ions 
contribute equally. The remaining contributions come from the 4/-states of Cel 
and Ce2, respectively (cf. Fig.|3]). The contributions of all cerium cations is 
consistent with the small tails of the 4/-states in the PDOS of Fig. WO^) at zero 
energy. In fact the entire vb contains contributions of 4/-states which decrease 
towards lower energies. 

The calculated GGA+U band structure yields insight into the optical properties of 
a-Ce2S3 and its color. The dispersion at the lower edge of the cb at k = F causes 
the rather small direct band gap of 0.8 eV. However, this gap is between states of 
Cel and Ce2. Since the color of the cerium sesquisulfides is determined by the 
excitations 4f^5d^ 4f^5d^ taking place on the same site we have to distinguish 
between two optical transitions. The orbital-projected band-structure of Fig. [3] 
reveals that the smallest Af-5d gaps of Cel and Ce2 are about 1.0 and 1.7 eV, 
respectively. Consequently, the reason why a-Ce2S3 has a brownish-black color is 
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twofold. The cerium cations Ce2 absorb photons with energies larger than 1 .7 eV 
which would lead to a dark red color. The smaller gap of ~ 1.0 eV (Cel cations) is 
responsible that all radiation of the visible range is absorbed. However, the 
absorption of radiation with an energy between 0.7 and 1.8 eV (near infrared and 
red light) is relatively low because the density of states is rather small in this 
region [compare Fig. [2tb)]. This effect leads to the black shade. Dark red mixed 
with black yields a brown-black hue. 



4.2 y-Phase 



Table |2] summarizes lattice parameters and selected interionic distances calculated 



for y-Ce2S3 and compares the values with that of the x-ray structure of Ref. Il34n . 
Again, the chosen value of U' yields a calculated unit-cell volume which is less 
than 1% larger than that determined from experiment. The deviation between the 
calculated and measured average Ce-S bond lengths is negligibly small. The 
smallest calculated Ce-Ce distance which is 2% smaller than that determined 
from the x-ray structure. The reason for the cell expansion is a growth of the voids 
located on the cation sublattice [cf. red spheres of Fig.[TJc)]. The calculated 
structure has distances between the center of the voids, Q, and the nearest sulfur 
anions which are between 3.1% and 4.0% larger than those of the experimental 
structure. 

The I42d structure of Fig. [U^c) has three symmetry inequivalent cerium cations 
and sulfur anions, respectively. Two of the symmetry inequivalent cerium cations 
have clearly distinguishable bond distances of r(Ce-S)~2.90 and r(Ce-S')~3.10 
A. However, the third kind of cations has an almost continuous spectrum of bond 
length to its eight nearest sulfur anions in the range of 2.90 and 3.30 A. Due to the 
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vacancies two of the three symmetry inequivalent sulfur anions are bound to five 
cerium cations only. All the penta-coordinated sulfur anions have three Ce-S 
bonds and two Ce-S' bonds [cf. Fig.fHb)]. The six-fold coordinated sulfur anions 
have Ce-S distances, which are between 1 and 2% larger than the calculated 
average bond lengths of Tab. |2l The calculated energy of formation for the y-phase 
is AEform = —814 kJ mol~^ Ce2S3. 

Figure |4] depicts the band structure and PDOS calculated for the fully optimized 
Cei6S24 supercell which is the primitive cell of that shown in Fig.fljc). 
Comparison with Fig.[2]reveals similarities between the electronic structures of 
the a-phase and the y-modification, but also significant differences which will be 
addressed in more details in Sec. \5\ Up to an energy of — 1 .0 eV, the vb is mainly 
formed by S 3p-orbitals with small contributions of Ce 4/- states. Between —1.2 
and —0.5 eV the vb has increasing admixtures of Ce 4/-states. The region of the 
occupied states above —0.5 eV (directly on top of the vb) is dominated by the 
LHB with significant components of S 3p- states. The difference to the PDOS of 
the a-phase is that the LHB of y-Cei S3 has one maximum only at around —0.4 
eV. The width of the LHB is around 0.5 eV without any distinct dispersion. The 
location of the 4/ LHB is in good agreement with available photoelectron 
spectroscopy data, but inferring its width from these experiments is difficult given 



the resolution [|l5|] . 



The single maximum of the LHB in the PDOS is evidence for cerium cations 
which are electronically and magnetically equivalent. This is consistent with the 
fact that the calculations yield for each cerium cation occupation numbers of 
Hf = 1.14. Occupation numbers larger than one indicate admixtures of other band 
states, mainly S Bp-orbitals. As in the case of a-Ce2S3 the cerium cations of the 
y-phase have a local effective magnetic moment of 2.6 /jb- 
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The cb has predominantly Ce 5d character with small admixtures of S 3jc-states 
and has its largest dispersion around k = P. At the same k-point is the smallest 
direct gap between the maximum of the LHB and the cb minimum of around 1.8 
eV which is in fact in the range of measured optical gaps between 1.8 and 2.0 eV 
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5211 . Thus the present calculations on Y-Ce2S3 result in a dark red color. 
The rather weak dispersion of both the LHB and of the Ce 5d bands at the bottom 
of the cb [see Fig. Ufa)] implies a rather sharp absorption edge. This is in contrast 
with the a-allotrope in which both the LHB and the 5d cb are more disperse. The 
smallest gap between the vb and the cb is 2.3 eV, also at k = P. The minimum of 
the lower edge of the UHB is located approximately 2.8 eV above zero energy. 



5 Comparison between the Ce2S3 modifications 

Both allotropes have crystal structures with a relatively high average coordination 
of eight sulfur atoms around each Ce atom. The mass density of the two Ce2S3 
phases calculated from the optimized crystal structures is similar: 5.10 Mg m^ 
(a-phase) and 5.13 Mg m^^ (y-phase). For comparison the experimentally 
determined mass densities for a-CeiSa and Y-Ce2S3 are 4.95 and 5.18 Mg m^^ 



ilOn . The calculations predict that the electronic energy of the a-phase is 6 
kJ mol^^ Ce2S3 lower (i.e. more stable) than that of the y-phase. This is in 
agreement with the fact that a-Ce2S3 is the low temperature configuration. In 
order to be more stable at higher temperatures the y-phase must have a larger 
entropy than the a-modification. This is plausible because the y-phase contains 
vacancies on the cationic sublattice. The voids might cause more low energy 
lattice phonons in the y-phase than in the a-phase. A lowering of phonon 
frequencies leads to a higher vibrational entropy. Furthermore, the distribution of 
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voids in the lattice adds a configurational entropy, which also tends to stabilize the 
y-phase at elevated temperatures. 

Figure [5] compares the partial density of states of a-Ce2S3 and y-CeiSa in the 
energy range of —37 to 10 eV. Overall, the electronic structures of the two 
allotropes have common features. The bonding between Ce and S atoms has 
significant covalent components. It is dominated by the interaction between Ce 
5d- and S 3/)-orbitals. The centers of the bonding and anti-bonding p-d bands are 
approximately 7 eV apart in both compounds (cf. Fig. [5]). 

In both modifications, the effective one-particle energy of this LHB is at the top of 
the S Bp vb. This leads to a hybridization between Ce 4/- and S 3p-orbitals. 
Electronic excitations in the visible part of the spectrum involve transitions from 
the occupied LHB into the cb, which is dominated by Ce BJ-levels. 

However, the present calculations reveal a major distinction between the a— and 
y— allotropes, due to the structural differences. The presence of two different 
coordination environments of the cerium cations in the a-phase causes a splitting 
near the minimum of the cb and also within the occupied 4/- states at the top of the 
vb (Figs. [2] and [5]). As a consequence, the optical transitions at 1.8 eV responsible 
for the red color are superimposed by transitions, which start around 0.8 eV. 
Furthermore, this distinction between the two types of Ce cations in the a-phase is 
noticeable also in the Ce 55- and Ce Sp-semi-core states at —34 and —17 eV, 
respectively (Fig. [5]). Both states are split by about 0.5 eV whereby the 5s and 
5/?-states of Cel are located at lower energies than that of Ce2. This implies that 
the effective potential for electrons around Cel is more attractive than around Ce2. 
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6 Summary and conclusions 



The present GGA+U calculations yield a very good description of the crystal 
structures and of the chromatic properties of a-Ce2S3 and Y-Ce2S3 involving 
4/-electrons explicitely. With a Coulomb repulsion parameter of U' =4 eV the 
calculated cell volume and most of the interionic distances for both Ce2S3 phases 
are less than 1% larger than that determined from x-ray diffraction measurements. 
This is comparable to the accuracy of standard GGA calculations as applied to 
systems without /-electrons. This is a major improvement over standard GGA 
calculations, which underestimate the cell volume of CeiSs derivatives by as 



much as 6% and fail to predict their semiconducting character OOll 



With the same U' parameter the present calculations yield effective one-electron 
energies, which allow a quantitative interpretation of optical transitions and thus 
accounts for the colors of the two Ce2S3 allotropes in a satisfactory manner. The 
vb and cb of both allotropes are dominated by S 3p and Ce 5<i-states, respectively. 
In both phases the occupied Ce 4/-states are energetically located between the vb 
and cb, whereby the smallest gap between the occupied Ce 4/- and the empty 
5J-states (A4y^_5rf) is direct and located at k = P. 

For the a-phase the calculations yield /^Af-5d = 0.8 eV and a rather small density 
of states in the range of 0.8-1.5 eV above the vb. The small density of states is due 
to three low-lying unoccupied bands with strong dispersion and corresponds to 
low absorption of radiation. The dispersion is related to the (Cel)S8 coordination 
polyhedron and it is the reason for the smaller optical band gap and explains the 
brownish-black hue a-Ce2S3. The dark red color (burgundy) of the Y-Ce2S3 is due 
to A4y^_5rf=1.8 eV and a rather steep increase of the density of states around the 
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lower edge of the cb. The y-allotrope is charcaterized by somewhat less dispersive 
bands and hence a sharper absorption edge. 

Thus, the present calculations provide an explanation for the difference of 1 eV in 
the optical gaps of a-Ce2S3 and y-CeiSs. The main factor is the presence of two 
types of Ce cations in the a-phase with two different coordination environments 
and the related difference of 3% between the average Ce-S bond distances in the 
two allotropes. Despite the different sulfur environments the calculations yield for 
cerium cations effective local magnetic moments of 2.6 hb which is in agreement 
with magnetic measuremts. 

It is gratifying that the computed electronic energy of the low-temperature 
modification (a-Ce2S3) is lower by 6 kJ mol^ (per formula unit) than that of the 
high-temperature modification (Y-Ce2S3). At low temperatures, the Gibbs free 
energy, which determines phase stabilities, is therefore lower for the a-phase. At 
elevated temperatures, the presence of vacancies in the y-phase is expected to 
enhance the vibrational and configurational entropy, thus stabilizing the 
high-temperature phase. 

The present GGA-i-U approach offers a satisfactory description of the Ce 4/-states 
in terms of structural and optical properties. However, the introduction of the 
on-site Coulomb repulsion does not cure the drawback of standard DFT methods 
of underestimating the gap between the itinerant S 3p- and the Ce SJ-band states. 
The smallest calculated direct gaps between the vb and the cb of y-Ce2S3 is 
computed to be 2.3 eV. This is approximately 1 eV smaller than the value 
determined from optical, photoconductivity, and photoluminescence spectra. The 
same band gap calculated for a-Ce2S3 is about 1.3 eV and probably also 
underestimated by 1 eV. A remedy of this problem might be the combination of 
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GGA+U and GW techniques ||53ll. 



In conclusion, the present work demonstrates that the GGA+U approach is 
successful in describing quantitatively the structural and electronic properties of 
a-Ce2S3 and Y-Ce2S3. The GGA+U approach overcomes a major shortcoming of 
standard GGA band structure calculations, which incorrectly yields a metallic 
bandstructure with 4/-states at the Fermi level. Instead, we found Ce2S3 to be an 
/-electron Mott insulator, in which strong correlation effects open up a Mott gap 
within the 4/-states, between a LHB and an UHB. The LHB is found to play a key 
role in the physical properties of these compunds, since the optical transitions 
responsible for their color are between the top of the LHB and the bottom of the 
5d cb. Furthermore, subtle differences in the calculated electronic structure of the 
two allotropes (most notably the presence of two inequivalent cerium sites in the 
a-modification, as well as more dispersive bands) give insight into the differences 
in the optical properties and color of these materials. The present GGA+U 
methodology applied here to RE compounds, could also be used to predict 
structural and opto-electronic properties of materials where RE elements are used 
as dopants. This capability will be particularly useful in the design and 
improvements of materials for lighting and display technologies with major 
benefits for our society and the environment. 
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List of Figures 



1 (color online) Structural features of cerium sesquisulfide 
allotropes. Blue and yellow spheres represent cerium cations and 
sulfur anions, respectively, (a) Coordination polyhedra of the 
symmetry inequivalent cerium cations Cel and Ce2 of a-Ce2S3. 
(b) Coordination polyhedron of Ce cations in Y-Ce2S3. The 
sulfur anions S and S' are within and out of the drawing plane, 
respectively, (c) Smallest conventional unit cell of Y-Ce2S3with 
a stoichiometric amount of cerium vacancies (red spheres). The 
composition of the super cell is Ce32S48 and has the highest 
symmetry possible (space group lAld). 

2 (color online) Electronic structure of a-Ce2S3 around the Fermi 
energy; (a) band structure and (b) partial DOS calculated for the 
fully optimized crystal structure. Arrows indicate relevant direct 
band gaps. The energy axes are labeled relative to the energy of the 
highest occupied band state. The latter is set to zero and is denoted 
by dashed lines. VB - valence band, CB - conduction band, UHB 
and LHB - upper and lower Hubbard bands, respectively. 

3 (color online) Contributions of atomic states per k-point to 
particular electronic bands of a-Ce2S3 between k = T and X as 
shown in Fig. |2l The thicker the electronic bands the larger the 
contribution of particular atomic states. Significant contributions 
of specific atomic states are indicated at k = P. Asterisks denote 
anti-bonding interactions. More details are given in the legend of 
Fig. [21 

4 (color online) Electronic structure of pure Y-Ce2S3 around the 
Fermi energy, (a) Band structure and (b) partial density of states 
calculated for the fully optimized system Cei6S24 (space group 
I Aid). More details are given in the legend of Fig. [2l 

5 (color online) Calculated partial density of states of two Ce2S3 
modifications in the energy range between —37 and 10 eV. 
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Fig. 1. (color online) Structural features of cerium sesquisulfide allotropes. Blue and yellow 
spheres represent cerium cations and sulfur anions, respectively, (a) Coordination polyhe- 
dra of the symmetry inequivalent cerium cations Cel and Ce2 of a-Ce2S3. (b) Coordination 
polyhedron of Ce cations in Y-Ce2S3. The sulfur anions S and S' are within and out of the 
drawing plane, respectively, (c) Smallest conventional unit cell of Y-Ce2S3with a stoichio- 
metric amount of cerium vacancies (red spheres). The composition of the super cell is 
Ce32S48 and has the highest symmetry possible (space group I42d). 
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Fig. 2. (color online) Electronic structure of a-Ce2S3 around the Fermi energy; (a) band 
structure and (b) partial DOS calculated for the fully optimized crystal structure. Arrows 
indicate relevant direct band gaps. The energy axes are labeled relative to the energy of the 
highest occupied band state. The latter is set to zero and is denoted by dashed lines. VB - 
valence band, CB - conduction band, UHB and LHB - upper and lower Hubbard bands, 
respectively. 
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Energy [eV] 

Fig. 3. (color online) Contributions of atomic states per k-point to particular electronic 
bands of a-Ce2S3 between k = T and X as shown in Fig. [2l The thicker the electronic 
bands the larger the contribution of particular atomic states. Significant contributions of 
specific atomic states ai^e indicated at k = P. Asterisks denote anti-bonding interactions. 
More details are given in the legend of Fig. [2l 
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Fig. 4. (color online) Electronic structure of pure Y-Ce2S3 around the Fermi energy, (a) 
Band structure and (b) partial density of states calculated for the fully optimized system 
Cei6S24 (space group lAld). More details are given in the legend of Fig.|2l 
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Fig. 5. (color online) Calculated partial density of states of two Ce2S3 modifications in the 
energy range between —37 and 10 eV. 
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List of Tables 



Calculated lattice parameters and selected interionic distances, 
Cel-X and Ce2-X, of a-Ce2S3 compared with data obtained from 



the x-ray diffraction measurements of Ref. nll\\ . The notation is 



according to that of Fig. [Ua). All distances are given in A. |34 

Comparison of calculated and measured lattice parameters and 
selected interionic distances of Y-Ce2S3. The calculated data have 
been obtained from the fully optimized supercell of Cei6S24 as 
depicted in Fig. [Itc). The experimental data has been taken from 



Ref. Il34ll . All values are given in A. |35 
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Table 1 

Calculated lattice parameters and selected interionic distances, Cel-X and Ce2-X, of a- 
CeaSs compared with data obtained from the x-ray diffraction measurements of Ref. lElll . 
The notation is according to that of Fig.[Tla). All distances are given in A. 

Lattice (space group Pnma) Interionic Distances 

X Cel Ce2 



Parameters 


Exp. 


Calc. 




Exp. 


Calc. 


Exp. 


Calc. 


a 


7.532 


7.565 


SI 






2.850 


2.872 


b 


4.097 


4.119 


SI' 


3.100, 


3.062, 


3.040 


3.056 










3.135 


3.161 






c 


15.728 


15.722 


S2 


2.983 


2.998 


2.883 


2.893 


AV/K,p [%] 




0.95 


S2' 


4.727 


4.761 


2.976 


2.988 








S3 


3.006 


3.018 


2.897 


2.893 








S3' 


2.850 


2.853 


4.696 


4.719 








Cel 


4.031 


4.029 


4.015 


4.035 








Ce2 


4.015 


4.035 


4.097 


4.119 
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Table 2 

Compaiison of calculated and measured lattice parameters and selected interionic distances 
of y-CeiSa- The calculated data have been obtained from the fully optimized supercell of 
Cei6S24 as depicted in Fig. [He). The experimental data has been taken from Ref. 13411 . All 
values are given in A. 



Lattice (Quasi cubic) 



Distances 



Parameter 


Exp. 


Calc. 




Exp. 


Calc. 


a 


8.637 


8.677 


''mm(Ce---Ce) 


4.040 


3.958 


c 


8.637 


8.641 «' 


(r(Ce-S)) ^) 


2.896 


2.898 


/Vexp X 100 




0.97 


(r(Ce-S')) 


3.087 


3.086 








KQ-S) ^) 


2.896 


2.988 








r(Q-S')^) 


3.087 


3.208 








rMm(Q---Q) ^' 


7.785 


7.779 



b) 



c/3 of the tetragonal supercell of Fig. [Tic). 

Q is the center of the voids on the cation sublattice. 



''^ The notation is according to Fig.[T}b). 
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